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SUMMARY 


Equations  governing  the  performance  of  charring  ablators^  subject  to  a 
diffusion- controlled  oxidation  mechanism  of  char  removal^  are  derived.  A  solu¬ 
tion  is  presented  for  the  case  of  a  const ant -enthalpy  level.  The  effects  of 
various  material  properties  and  environmental  parameters  are  discussed.  A'  ■ 


INTRODUCTION 


The  charring  plastics  provide  a  practical  method  for  protecting  the  inte¬ 
rior  of  a  reentry  vehicle  from  the  aerodynamic  heating  encountered  during 
reentry.  Much  effort  has  been  directed  toward  understanding  the  performance  of 
these  materials.  Analytical  models  for  charring  ablators  are  presented  in  ref¬ 
erences  1  and  2.  These  references  also  present  estimates  of  the  weight  of  the 
charring  ablator  which  is  required  to  provide  adequate  thermal  protection  as  a 
function  of  heating  conditions  and  material  properties.  A  more  rigorous  ana- 
lyticaJ-  model  of  the  charring  ablator  is  presented  in  reference  3*  Results  of 
some  typicad  calculations  are  given  in  references  3  to  5,* 

A  primary  problem  encountered  in  analyses  of  the  performance  of  charring 
ablators  is  the  removal  of  char  from  the  surface.  Aerodynamic  shear^  theimal 
stress^  differential  pressure^  and  oxidation  have  been  suggested  as  mechanisms 
contributing  to  char  removal.  The  removal  of  char  as  a  result  of  pressure 
within  the  char^  caused  by  confinement  of  the  gaseous  products  of  pyrolysis^  is 
analyzed  in  reference  6.  The  effects  of  aerodynamic  shear  and  thermal  stress 
are  not  known  in  quantitative  terms.  Results  of  an  extensive  experimental 
investigation  of  the  effects  of  oxidation  on  char  removal  and  surface  heating 
are  presented  in  reference  7*  the  test  environment  employed  in  the  inves¬ 

tigation  (ref.  7)^  oxidation  appeared  to  be  the  only  mechanism  involved  in  char 
removal.  In  the  light  of  the  results  presented  in  reference  7^  it  appears  that 
oxidation  effects  should  be  incorporated  in  the  analysis  of  charring  ablators. 

The  differential  equation  governing  the  performance  of  a  charring  ablator 
subject  to  a  diffusion  controlled  mechanism  of  oxidation  has  been  derived 
herein  and  is  solved  for  a  particular  case.  The  effects  of  material  properties 
and  environmental  factors  on  performance  are  discussed. 


SYMBOLS 


^  -v.-r0.6 

a  =  CeA^Le 

1  +  Tl(a^  +  -  H,) 

Aheff(l  +  ari) 

C  concentration  of  oxygen 

Ce  concentration  of  oxygen  external  to  boundary  layer 

f  volatile  fraction  in  un charred  material 

Z^h  heat  of  combustion  per  unit  weight  of  char  consumed 

^eff  effective  heat  of  pyrolysis 

Z^p  heat  of  pyrolysis 

H  total  enthalpy 

k  thermal  conductivity 

m  rate  of  mass  transfer 

m^  rate  of  char  removal 


N3^  Schmidt  number 

q.p  heating  of  surface  resulting  from  combustion 


hot- wall  convective  heating  rate 
1(3^0  hot-wall  convective  heating  rate  with  no  mass  transfer 

radiant  heating  rate 
Q  heat  input 
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total  convective  heat  input 
total  radiant  heat  input 


absolute  temperature 

average  surface  temperature  during  reentry- 
boundary- layer  velocity  parallel  to  surface 
boundary- layer  velocity  normal  to  surface 

thermal- protect ion  weight  per  unit  area  (weight  of  material  degraded) 
char  weight 

boundary- layer  coordinate  parallel  to  surface 
boundary- layer  coordinate  normal  to  surface 


v'^ref  ^p  )('^ref  "^p 


char  thicknes 


emissivity 

transpiration  factor 


weight  of  char  removed  per  unit  weight  of  oxygen  diffusing  to  surface 

viscosity 

char  density 

density  of  boundary- layer  fluid 


Stefan-Boltzmann  constant 


Subscripts: 


external  to  boundary 
initial  value 


extremum  value 


max 


maximum  value 


0  value  in  absence  of  mass  transfer 

p  pyrolysis 

r  radiant 

ref  reference  value 

s  surface 

w  at  the  vail 


AMLYSIS 


(Reradiatian) 

f 


j  “‘c 


Char  lays  I' 

?  i  t 

1 

Uti char  red 
material 

Front 

surface 


Pyrolysis 

interface 


An  analytical  model  of  the  charring  ablator  is  presented  in  reference  1. 

This  model  is  shown  schematically  in  figure  1.  The  outer  surface  is  subjected 

to  aerodynamic  heating.  The  char  pro¬ 
vides  an  insulating  layer^  and  the  heat 
passing  through  this  layer  is  absorbed 
by  pyrolysis  at  the  interface  between  the 
char  layer  and  the  uncharred  material. 

The  surface  of  the  char  re-radiates  a 
significant  fraction  of  the  heat  input. 
The  gases  generated  by  pyrolysis  tran¬ 
spire  through  the  char  layer  and  are 
injected  into  the  boundary  layer.  The 
gases  are  heated  as  they  pass  through 
the  char  layer,  and  this  removal  of  heat 
from  the  char  layer  reduces  the  quantity 
of  heat  conducted  to  the  pyrolysis  inter¬ 
face.  When  these  gases  are  injected  into 
the  boundary  layer,  the  convective  heat  transfer  to  the  surface  is  reduced. 

This  reduction  in  convective  heat  as  a  result  of  mass  injection  is  the  same 
effect  as  that  obtained  with  simple  subliming  ablators.  In  addition  to  the 
gases  produced  by  pyrolysis,  a  residue  remains  at  the  interface  and  adds  to  the 
thickness  of  the  char  layer.  This  process  tends  to  Increase  the  thickness  of 
the  char  layer.  At  the  same  time,  char  removal  may  occur  as  a  result  of  thermal 
chemical,  or  mechanical  processes,  and  tends  to  decrease  the  thickness  of  the 
char  layer. 


Figure  1.-  Schematic  diagram  of  charring 
ablator. 


thermal  behavior  of  charring  ablators  is  very  complex  and^  in  general 
niimerical  methods  are  required  to  solve  the  governing  equations.  However^  if 
some  simplifying  assumptions  are  made^  approximate  expressions  can  be  obtained 
for  the  char  accumulation  and  the  pyrolysis  rate.  Such  approximate  results  are 
compared  with  the  numerical^ results  of  reference  1,  and  the  agreement  is  satis- 
f actox^  .  The  effects  of  oxidation  on  char  removal  and  heat  input  can  be 
included  i.n  analyses  similar  to  that  of  reference  1  without  any  particular  dif¬ 
ficulty.  However^  equations  must  be  obtained  for  the  rate  at  which  oxidation 
removes  char  and  the  rate  at  which  oxidation  adds  heat  to  the  surface. 


Oxidation 


It  is  sliown  in  roiGroncc  8  that  ttiG  rate  of  oxidation  of  porous  carbon  in 
air  is  a  diffusion  controlled  process,  even  at  relatively  low  temperatures. 
Therefore,  the  rate  of  oxidation  is  controlled  by  the  rate  at  which  oxygen  dif¬ 
fuses  from  the  external  stream  to  the  surface.  A  detailed  study  of  the  removal 

of  char  by  oxidation  is  presented  in  reference  7 •  that  study  it  is  concluded 

that  char  removal  occurs  primarily  as  a  result  of  oxidation.  With  sufficiently 
surface  temperatures  (and  hence  high  reaction  rates),  the  oxygen  diffusing 
to  the  surface  reacts  with  the  char  rather  than  with  the  products  of  pyrolysis. 
The  msiximum  rate  at  which  oxygen  diffuses  to  the  surface  is 


m(02) 

The  derivation  of  equation  (l)  is 

If  A  is  the  weight  of  char 
the  rate  of  char  removal  is 


C 


(1) 


He  - 

outlined  in  the  appendix. 

removed  per  unit  weight  of  available  oxygen^ 


mn  = 


-  H 


w 


^  ^.0.6- 
^e^Le  ^ 


^c 
He  “ 


(2) 


If  Ah  is  the  quantity  of  heat  added  to  the  surface  per  unit  weight  of  char 
consumed^  the  rate  at  which  the  surface  is  heated  by  combustion  is 


^c 
He  - 


a  Ah 


(3) 


The  evaluation  of  A  and  Ah  is  a  complicated  aero- thermal- chemical 
problem^  the  solution  of  which  is  beyond  the  scope  of  the  present  paper.  How¬ 
ever^  limiting  values  can  be  indicated.  The  maximum  value  of  A  is  0-75^ 
which  corresponds  to  the  ratio  of  the  molecular  weights  of  carbon  and  oxygen 
when  the  product  of  oxidation  is  carbon  monoxide.  For  this  reaction^  the  value 
of  Ah  is  approximately  4^000  Btu/lb.  Carbon  monoxide  is  the  most  probable 
product  of  a  reaction  occurring  in  an  oxygen-deficient  environment,  and  the 
formation  of  it  is  consistent  with  the  results  of  reference  7* 

The  model  of  a  charring  ablator  described  previously  is  adequate  for  numer¬ 
ical  analysis.  With  certain  simplifications,  the  model  can  be  reduced  to  a  form 
for  which  analytical  solutions  can  be  obtained.  First,  it  is  assumed  that  the 
specific  heat  of  the  char  layer  is  negligible.  Solutions  based  on  this  assump¬ 
tion  are  compared  with  numerical  solutions  in  reference  1  and  the  approximate 
solution  provides  an  adequate  representation  of  phenomena  within  the  char, 
including  char  formation.  A  second  assumption  is  that  the  effective  conduct iv-r 
ity  of  the  char  layer  is  a  function  of  temperature  as  follows: 
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Because  the  specific  heat  of  the  char  layer  is  negligible^  the  heat  transferred 
across  the  char  layer  (from  the  surface  to  the  reaction  zone)  is  equal  to 
k(Ts^Tp) 

It  is  expected  that  the  thermal  conductivity  of  a  porous  material 
^c 

such  as  a  char  layer  "will  increase  with  increasing  temperatiire  because  of  the 
greater  importance  of  radiation  across  the  void  spaces.  The  surface 
re-radiates  energy  at  a  rate  proportional  to  its  fourth  power;  therefore^  rela¬ 
tively  large  changes  in  the  heat  input  are  required  to  obtain  significant 
changes  in  the  surface  temperature.  Therefore^  by  properly  selecting  the  refer¬ 
ence  temperature  T^^P^  the  error  introduced  by  using  equation  can  be  mini¬ 
mized.  The  actual  error  cannot  be  estimated  because  no  data  are  presently 
available . 

The  model  of  a  charring  ablator  presented  in  reference  1  can  be  readily 
extended  to  the  case  in  which  char  removal  and  surface  heating  occur  as  a  result 
of  oxidation.  The  unknown  quantities  are  the  surface  temperature^  the  rates  of 
pyrolysis  and  char  removal^  the  convective  and  combustive  heating  rates^  and 
the  weight  of  char  accumulated.  From  these  quantities^  the  weight  of  material 
degraded  can  be  derived. 

The  approach  which  will  be  used  is  as  follows .  Energy  balances  at  the 
heated  surface  and  at  the  pyrolysis  interface  are  used  to  eliminate  the  surface 
temperature.  An  approximate  boundary- layer  solution  is  used  to  express  the  con¬ 
vective  heating  rate  in  terms  of  the  convective  heating  rate  with  no  mass  trans¬ 
fer  and  the  mass-transfer  rates.  The  char- removal  rate  and  the  combustive- 
heating  rate  are  related  to  the  convective-heating  rate  by  equations  (2)  and  (j). 
An  equation  for  the  weight  of  char  accvunulated  is  obtained  by  noting  the  com¬ 
bined  effects  of  char  formation  and  char  removal.  A  differential  equation 
relating  the  weight  of  char  accumulated  and  the  convective  heat  input  with  no 
mass  transfer  is  then  derived.  This  equation,  together  with  the  relation 
between  char  weight  and  amount  of  material  pyrolyzed,  is  used  to  relate  the 
weight  of  material  degraded  and  the  convective-heat  input  with  no  mass  transfer. 

A  heat  balance  at  the  outer  surface  yields  the  following  equation: 

qc  +  q^  +  qj.  =  oeTs^  +  |^(Ts  -  Tpj 

Convective  Combustion  Radiant  Re-radiation  Heat  conducted 
heating  heating  heating  from  surface 

The  heat  conducted  from  the  surface  is  absorbed  by  pyrolysis  and  by  heating 
the  gaseous  products  of  pyrolysis  to  the  surface  temperature.  The  rate  of 
pyrolysis  is 
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(6) 


-  'tp) 


where 


Q 

Ahgff  =  +  "^('^s  “  '^p) 

is  an  effective  heat  of  pyrolysis  which  includes  the  heat  absorbed  by  the  gases 
as  they  transpire  through  the  char.  Equation  (7)  is  derived  in  reference  2. 
Equations  (5)^  (^)  combined  to  eliminate  T 


^r  +  ^c  =  ^ 


(8) 


where 


p  ^  _ frefP _ 

(^ref^  '^p^)  ('^ref  '*■  '^p) 


and  it  is  assumed  that 


Tp^  «  Tg^ 


The  convective-heating  rate  can  be  obtained  from  the  relation 

.  \®e  “  ^ 


‘Ic  =  ^Ic.o 


-  n(mc  +  “p)”Ti 


c,0 


(9) 


where  q^  q  is  the  hot-wall  convective -heating  rate  with  no  mass  transfer. 

Equation  (9)  is  obtained  by  fitting  a  straight  line  to  the  boundary- layer  solu¬ 
tions  given  in  reference  9-  It  is  valid  only  when  the  term 


V Hg  - 

'  ^c,o~ 


This  condition  is  almost  always  satisfied  when  all  the  heat  input  is  convective. 
A  higher  order  approximation  may  be  required  if  high  radiant-heat  inputs  are 
also  incident  on  the  surface.  It  is  assumed  in  equation  (9)  that  the  char  which 
is  oxidized  and  the  pyrolysis  products  which  are  injected  into  the  boundary 
layer  are  equally  effective  in  blocking  the  convective  heat  input;  this  assump¬ 
tion  explains  the  multiplication  of  me  Q^nd  nip  by  the  same  value  of  t]. 

Actually,  the  blocking  effectiveness  depends  on  the  molecular  weight  of  the 
injected  gases  and  probably  differs  for  the  products  of  combustion  and  pyrol¬ 
ysis.  However,  to  the  present  order  of  approximation  this  difference  may  be 
neglected. 
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From  equations  (2)  and  (9),  the  rate  of  oxidation  of  the  char 


1  +  aTi\He  -  H 


and  similarly  the  heat  transfer  to  the  surface  (in  terms  of  the  heat  transfer 
with  no  mass  transfer  and  the  rate  of  pyrolysis)  is 

q  _  (Sg  -  Hy)mg  _  *^0,0  ~  l^l^e  ~  ^w)  . 

^  1  +  arj 

From  ec^uations  (5)  and  (ll)^  the  rate  at  which  the  surface  is  heated  as  a 
result  of  combustion  is 


®e  “  ^w 


1  +  arj 


To  this  pointy  a  sufficient  number  of  relations  has  been  developed  to 
express  all  the  unknown  variables  in  terms  of  the  rate  of  pyrolysis^  the  weight 
of  char  accumulated^  and  the  heating  rate  with  no  mass  transfer.  From  consider- 
ation  of  the  geometric  effects  of  pyrolysis  and  char  removal^  the  weight  of  char 
accumulated  is  related  to  the  rates  of  pyrolysis  and  char  removal  as  follows: 

1  -  f  . 

iT  =  “7“  “P  ■ 

which  states  that  the  rate  of  change  of  char  weight  is  equal  to  the  difference 
between  the  rate  at  which  char  is  formed  and  the  rate  at  which  char  is  removed. 
Equation  (lO)  can  be  used  to  eliminate  m^  from  equation  (15).  Thus^ 

dWc  .  fl  -  f  +  an1 


f(l  +  aq) J  (1  +  aq)(He  -  H^) 


Equations  (8),  (ll),  (12),  and  (llj-)  can  be  combined  to  yield  the  fol¬ 
lowing  equation: 


^10,0  a  Ah  +  Eg- 


^e  *"  ^w  ^eff(^  +  arj)  Ah^^^ 


1  -  f  +  ai,[_  dt  H,.  -  hJ  U  i5h.f,(l  +  an) 


The  quantities  general^  known  functions 

of  time.  In  an  important  special  case  to  be  discussed  later^  that  is^ 


laboratory  evaluation  of  charring  ablators^  is  a  function  of  time  whereas 

Ic  o^  ^r^  He  -  are  constants.  Equation  (I5)  is  a  nonlinear  first- 

order  ordinary  differential  equation  which  can  be  solved  by  numerical  methods. 
When  the  weight  of  char  has  been  determined^  equation  (l4)  can  be  integrated  to 
determine  the  mass  of  the  gaseous  products  formed  by  pyrolysis.  The  weight  of 
material  degraded  to  yield  gaseous  products  and  char  is 

w  =  ^  +  Wc^i  (16) 


where  is  the  char  initial  weight. 

The  equation  governing  the  performance  of  a  charring  ablator  (eq.  (I5)) 
cannot  be  integrated  in  general  terras;  however^  certain  significant  information 
can  be  deduced  from  the  equa¬ 
tion  without  obtaining  a  solu¬ 
tion.  In  addition^  a  solution 
can  be  obtained  for  a  special 
case. 

Extreranm  Char  Weight 

When  the  derivative  of  the  W^ 

char  weight  with  respect  to 
time  vanishes^  the  char  weight 
is  an  extremum.  It  will  be  seen 
later  that  the  char  weight  is  a 
raonotonic  function  of  heat  input 
if  other  conditions  are  constant. 

Therefore^  the  extremum  may  be 
either  a  maximum  or  a  minimum 
depending  on  whether  the  char 
weight  is  lesser  or  greater  than 
the  extremum  char  weight,  as 
shown  in  sketch  (a). 

From  equation  (I5),  the  extremum  char  weight  is  found  to  be 

-  1  (IT) 


It  should  be  noted  that,  in  general,  is  a  function  of  time. 


w\l  -  f 


Ahgff  a  Ahgff 


(Hg  -  (1  -  f  +  ari) 


Hie  extremum  char  weight  is  an  iii5)ortant  combination  of  material  properties 
and  environmental  parameters,  entirely  aside  from  its  interpretation  as  a  lim¬ 
iting  value.  The  extremum  char  weight  depends  on  the  heating  only  through  the 
ratio  of  radiative  to  convective  heating  rates  ^r/lg  o-’  increases  as  the 

value  of  this  ratio  increases.  The  volatile  fraction  and  the  stream  enthalpy 
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strongly  influenced  the  extremum  char  weight,  as  shown  in  figures  2  and  3.  The 
effects  of  oxygen  concentration  and  radiant  heating  are  discussed  later. 


Figure  2.-  Effect  of  volatile  fraction  on  extremxim 
char  weight.  A  =  O.75;  Cq  =  0.20;  p  =  0.10. 


When  the  char  weight  approaches 
the  extreminn  char  weighty  dW^/dt 
approaches  zero^  from  the  defini¬ 
tion  of  extremum  char  weight.  When 
dWc 

=  Oj  the  rate  of  pyrolysis  is^ 
from  equation  (l4)^ 


1  -  f  +  ar]  He  - 


For  this  condition^  the  rate  of 
pyrolysis  is  independent  of 

and  The  only  material 

property  affecting  the  pyrolysis 
rate  is  the  volatile  fraction  f. 
This  result  can  be  understood  when 

it  IS  noted  that,  with  — =  0, 

dt  ^ 

the  rate  at  which  char  is  formed 
(which  is  directly  proportional  to 
the  rate  of  pyrolysis)  is  equal  to 
the  rate  of  char  removal.  However^ 
the  rate  at  which  char  is  removed 
depends  only  on  the  rate  at  which 
oxygen  diffuses  to  the  surface.  It 


should  be  noted  that  the  time 

required  to  approach  the  extremum  char  thickness  increases  as  the  extremum  char 
thickness  increases.  The  extremum  char  thickness  is  a  function  of  the  material 
properties  as  well  as  of  the  environmental  conditions;  therefore^  the  extent  to 
which  a  given  test  will  be  influenced  by  the  char  weight  approaching  its  extremum 
value  is  related  to  the  properties  of  the  material. 


Caution  must  be  exercised  in  the  interpretation  of  test  results  when  the 
char-formation  characteristics  in  the  test  environment  do  not  simulate  those 
that  exist  in  the  flight  environment.  When  the  char  thickness  approaches  its 
extremum  value  (for  example^  in  low  enthalpy  tests the  response  of  the  material 
is  governed  by  relations  which  are  quite  different  from  those  that  apply  in 
general.  The  only  material  characteristic  which  affects  material  performance 
under  such  conditions  is  the  volatile  fraction. 


Solution  for  Constant  Enthalpy 

The  most  important  single  heating  condition  is  that  for  which  the  stream 
enthalpy  is  constant  and  for  this  condition  equation  (I5)  can  be  integrated  in 
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closed  form,  provided  q.r/lc,o  *^e  constant.  All  of  these  conditions 

are  typical  of  tests  conducted  in  ground  test  facilities. 

Equation  (15)  can  be  simplified 

to  the  following  form  by  using  equa-  _ _ _ _ 

tion  (17):  “  '  0^2  / 

a  ^c,m  ^c\  1  / 

dW„  1  +  aql  ^  P/Hg  -  8  / _ o.if _ 


dWg  1  +  aq  \  p  Py  “  ^w 

d^  "  Wc  ,  ^  ^  +  He  - 

p  ^eff(l  ®''n) 


If  the  stream  enthalpy  is  constant  and 
any  radiant  heating  is  proportional  to 
the  convective  heating,  equation  (19) 
is  a  linear  first-order  differential 
equation  in  Qq  vlth  as  the  inde¬ 

pendent  variable.  The  solution  of  this 
equation  is 


-  1  + 


f  1 
0.2  / 

n.U  . 

/ 

/ 

/ 

// 

a 

1  +  ari  (He  -  H-w)Wc^in 


where 


A  =  1  + 


(20a) 

Ti(aAh  +  Hg  -H^) 
Aheff(l  +  s-'h) 


Figure  3.-  Effect  of  enthalpy  on  extremum 
char  weight,  =  0. 


(20b) 


The  amount  of  pyrolysis  that  has  occurred  can  be  determined  as  a  function  of  the 
heat  input  and  the  char  thickness  by  integrating  equation  (ik) 

a 

Wg  -  Wg  ,•  +  - - 

^  ^  1  +  at]  H0  -  (21) 

“  1  -  f  +  aq 

f(l  +  aq) 

The  weight  of  material  which  has  been  thermally  degraded  (that  is,  pyro- 
lyzed  to  yield  gaseous  products  or  char)  is,  from  equations  (l6)  and  (21), 

a  *^c 

^c  *"  ^c  i 

’  1  +  aq  Hg  -  ^  ^ 

W  =  - - - - -  +  Wg  i  (22) 

1  -  f  +  aT] 

1  +  aq 
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Because  of  the  transcendental  character  of  equation  (20)^  the  required  weight 
of  material  cannot  be  expressed  explicitly  as  a  function  of  the  heat  input. 
However^  equation  (20)  can  be  used  to  determine  Wq  as  a  function  of  Q^. 

These  results  can  be  substituted  into  equation  (22)  to  calculate  the  required 
weight . 


Solution  for  No  Char  Removal 


For  the  case  in  which  no  char  removal  occurs  (a  =  O)^  equation  (15)  has 
the  following  form: 


dWc 

^  'n(®e  “  ^w) 

^Cj  0  ^  ^r 

dt 

P  Aheff 

^eff 

(23) 


This  form  of  the  equation  would  apply  to  a  material  having  characteristics 
similar  to  a  charring  ablator  except  that  the  surface  does  not  oxidize  and  is 
not  removed  by  mechanical  forces.  An  impregnated  ceramic  shield  might  provide 
an  example  of  this  type  of  performance.  This  equation  would  also  apply  to 
reentry  into  a  planetary  atmosphere  containing  no  oxygen  or  other  chemical 
species  which  react  with  the  char. 


The  solution  of  equation  (23)  is 


^c,±  ^  ^  ^  ~  \)~|  2(1  -  f)  (^C  +  ^r. 

_  P  ■^eff  —  f  P  ^eff 


1/2 


W, 


Hli  +  1 


E  -  E  ) 
e  w/ 


Ah 


eff 


(2ifa) 


and  the  total  weight  of  degraded  material  is 


W  =  W 


c.i 


(24b) 


If  an  average  value  of  the  enthalpy  during  reentry  is  assumed,  this  result  is 
identical  to  that  presented  in  reference  1. 


DISCUSSION 


Certain  characteristic  features  of  char  formation  and  removal  have  been 
noted  previously.  In  general^  however^  the  equations  governing  the  performance 
of  charring  ablators  are  too  complex  to  yield  significant  information  from  a 
cursory  examination.  The  effects  of  material  properties  and  environmental 
parameters  on  system  performance  can  best  be  evaluated  on  the  basis  of  a  para¬ 
metric  study.  All  the  resxlLts  presented  in  figures  4  to  9  are  based  on  heating 
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in  an  airstream.  Therefore,  with  a  given  heat  input  but  with  different^  enthal¬ 
pies  the  difference  between  the  two  results  indicates  the  extent  to  which  a 
test  at  one  condition  fails  to  simulate  performance  at  the  other  condition;  that 
is,  it  indicates  the  difference  which  results  from  failure  to  consider  the  varia- 
tions  in  the  two  test  conditions. 

The  relation  between  w/p  and  — -§ -  is  independent  of  the  value 

P  Aheff 

of  p.  However,  the  required  weight  itself  is  not  independent  of  p.  For 
example,  if  p  and  Q  both  increase  by  a  factor  of  2,  the  required  weight  also 
increases  by  a  factor  of  2.  Little  information  is  available  concerning  the 
appropriate  value  of  p.  Based  on  a  comparison  of  calculated  res\ilts  and  a 
limited  number  of  test  results,  the  value  of  p  for  a  5O-5O  mixture  of  phenolic 
and  nylon  was  found  to  be  about  0.10. 


Thermal-Protection  Weight 


The  conditions  for  which  equation  (I9)  has  been  integrated  are  not  repre¬ 
sentative  of  an  actual  spacecraft  reentry.  During  such  a  reentry,  the  enthalpy 
decreases  from  its  Initial  value  to  a  negligible  value;  however,  a  significant 
fraction  of  the  total  heat  input  may  be  experienced  before  the  decrease  in 
enthalpy  is  appreciable.  In  tests  conducted  in  ground  facilities,  a  constant- 
enthalpy  level  is  normally  maintained  and  the  present  solution  should  be  valid. 
Therefore,  weights  calculated  from  equations  (20)  and  (22)  are  usefi^  in  eval¬ 
uating  the  effects  of  environmental  parameters  and  material  properties  on 
thermal-protection  requirements  and  also  in  determining  the  significance  of 
ground  test  results  as  they  apply  to  a  flight  environment. 

The  thermal-protection  weight  (eq.  (22))  is  plotted  as  a  function  of  the 
total  heat  input  in  figure  k.  Although  this  weight  does  not  include  any  allow¬ 
ance  for  insulation  or  for  factors  of  safety,  it  will  be  referred  to  as  the 
thermal-protection  weight.  The  required  weight  is  higher  at  the  low  enthal^^ 
level  than  at  the  high  enthalpy  level  for  all  values  of  total  heat  input.  This 
condition  results  from  the  greater  effect  of  oxidation  and  the  lesser  effect  of 
blocking  at  the  lower  enthalpy. 


Over  the  most  important  range  of  heat  input,  the  optimum  value  of 
between  0.4  and  0.6,  the  lower  value  of  f  being  associated  with  the  higher 
heat  input.  However,  for  values  of  f  between  0.4  and  0.6,  the  required  weight 
is  not  extremely  sensitive  to  the  particular  value  that  is  used. 


The  effect  of  total  heat  input  on  required  weight  is  shown  in  figure  5  for 
a  value  f  =  0.4.  With  this  linear  total  heat  input  scale,  the  decrease  in  the 
slope  of  the  weight  curve  with  increasing  total  heat  input  is  readily  observed. 
With  no  radiant  heat  input,  the  required  weight  is  mOre  strongly  influenced  by 
the  enthalpy  level.  If  one— half  of  the  total  heat  input  is  radiant 
(ir/^  o  “  ^)-»  required  weight  at  a  constant -enthalpy  level  is  less  than 

that  for  similar  total  heating  in  the  absence  of  radiation.  This  effect  is  much 
more  pronounced  at  the  low-enthalpy  level.  Available  test  results,  which  are 
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Figure  1)^.-  Effect  of  heat  input  on  required  weight  of  ablating  material. 


reported  on  the  hasis  of  temperature  rise  rather  than  on  the  weight  of  material 
degraded,  indicate  that  the  effect  of  radiant  heating  on  material  performance 


Figure  5.-  Effect  of  total  heating  on  required 
weight  of  ablating  material,  f  =  0.4. 


is  similar  to  an  equal  amount  of  con¬ 
vective  heating. 

The  effect  of  the  enthalpy  level  on 
required  weight  is  shown  explicitly  in 
figure  6  for  selected  values  of  total 
heat  input.  At  the  lower  enthalpies^ 
the  required  weight  is  very  sensitive  to 
the  enthalpy  value  whereas  at  higher 
enthalpies  the  required  weight  is  rela¬ 
tively  insensitive  to  the  enthalpy. 


Thermal-Protection  Effectiveness 

It  is  difficult  to  define  a  simple 
parameter  of  considerable  generality 
which  provides  a  measure  of  the  effec¬ 
tiveness  of  char-forming  materials  in 
the  sense  that  the  heat  of  ablation  pro¬ 
vides  such  a  measure  for  subliming 


ablators.  A  parameter  consisting 
of  a  total  heat  input  divided  by 
the  thermal-protection  weight 
required  to  accommodate  this  heat 
input,  which  has  been  used  in  a 
number  of  investigations,  appears 
to  provide  the  best  available 
characterization  of  charring 
ablators: 

Qc 

E  =  +  Qr 

The  weight  in  this  eq^uation  is 
that  required  to  assure  the  back 
surface  of  a  shield  subjected  to 
the  total  heat  input  ^r 

will  not  exceed  the  temperature 
of  pyrolysis.  Some  results  of 
the  present  study  are  presented 
in  terms  of  the  ratio  of  this 
parameter  to  figure  7  • 


The  effectiveness  increases  with  increasing  stream  enthalpy  and  is  greater 
for  high  total  heat  input  than  for  lower  heat  input.  The  values  shown  in  fig¬ 
ure  7  are  for  a  value  of  f  of  O.ij-  which  is  less  than  the  optimum  value  at  high 
enthalpies.  Therefore  the  estimated  system  effectiveness  in  figure  7  is  less 
than  the  maximum  theoretical  value.  At  higher  enthalpies  the  effectiveness  of 
charring  ablators  is  8  to  li-  times 
greater  than  the  effective  heat  of  i6 

pyrolysis  for  the  range  of  impor¬ 
tant  heat  loads. 

12 

Accumulation  of  Char 

E 

The  extremum  char  weight  is  a  8 

significant  parameter  in  the  per¬ 
formance  of  charring  ablators; 
however^  it  is  a  limiting  value 
and  not  necessarily  a  value  that 
will  be  attained  or  even  approached 
during  a  given  heating  cycle.  For 
example^  with  a  low  volatile  frac-  ^ 

tion  and  a  high  enthalpy  levels 
extremum  char  weights  in  excess  of 
20  pounds  per  square  foot  (assuming  Figure  T*-  Effect  of  enthalpy  level  on  effective- 

P  =  0.10)  are  predicted  from  fig-  ^^lating  material,  f  =  0.4;  <ir  =  0- 

ure  2.  This  value  is  greater  than 
the  estimated  total  thermal 
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protection  weight  requirements  for  the  stagnation  point  of  a  vehicle  entering 
the  earth’s  atmosphere  at  hyperbolic  velocity.  (See  ref.  4.)  Therefore^  it 
is  necessary  to  examine  actual  char  weights  in  relation  to  extremum  char  weight 

as  a  function  of  the  total 


heat  input.  The  char  weight 
as  a  fraction  of  the  extremum 
char  weight  is  plotted  in 
figure  8  as  a  function  of 
total  heat  inputs  with  zero 
initial  char  weight. 

At  the  lower  enthalpy 
level  (fig.  8(a) 

approaches  at  rela- 

tively  lov  total  heat  inputs^ 
particularly  when  the  vola¬ 
tile  fraction  is  high.  At 
the  higher  enthalpy  level 
(fig.  8(b) the  char  thick- 
ness  is  considerably  less 
than  its  maximum  value ^  even 
with  high  total  heat  inputs 


(a)  Low  enthalpy.  - 

Aheff 


and  a  high  volatile  fraction. 
Therefore^  at  low  enthalpy 
levels  the  performance  of  a 


charring  ablator  is  strongly 


H  —  H 

(b)  High  enthalpy.  — ^ - —  =  lO. 

^eff 

Figure  8.-  Effect  of  heat  input  on  the  relative  magni¬ 
tudes  of  char  weight  and  extremum  char  weight. 


influenced  by  the  char  thick¬ 
ness  approaching  its  maximum 
value ^  whereas  at  high 
enthalpy  levels^  the  approach 
of  the  char  weight  to  a  max¬ 
imum  value  is  barely  evident . 
When  radiant  heating  is  pres¬ 
ent^  the  char  weight  as  a 
fraction  of  its  maximimi  value 
is  greatly  reduced  for  given 
total  heat  input. 

That  ^cy^^c^max  ^ 
dominant  factor  in  the  per¬ 
formance  of  charring  ablators 
is  evident  from  the  basic 
equation  (eq.  (20)).  The 
total  heat  input  required  to 
produce  a  char  weight  (with 
Wc^i  =  0)^  that  is^  within 
l/e  of  the  maximum  char 
weight  is  plotted  in  fig¬ 
ure  9  as  a  function  of  f. 


An  increase  in  the  enthalpy 
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level  by  a  factor  of  5  results  in 
an  increase  by  a  factor  of  20  in 
the  heat  input  required  to  achieve 
the  given  value  of  ^c/^OymeiX- 


LIMITATIONS  AND  EXTENSIONS 
OF  ANALYSIS 


The  method  of  analysis  pre¬ 
sented  herein  is  based  on  a  dif¬ 
fusion  controlled  oxidation  of  the 
char  surface.  This  mechanism  has 
been  verified  in  a  series  of  tests 
of  50-percent -phenolic — 50-percent- 
nylon  specimens.  These  tests  were 
conducted  in  an  atmospheric- 
pressure  arc-heated  jet  with  a 
heating  rate  of  about 
180  Btu/ft^-sec.  For  tests  of 
other  materials  or  for  tests  in  a 
different  environment^  other 
mechanisms  of  char  removal  may  be 
important.  Even  when  oxidation 
is  the  only  mechanism  involved  in 
char  removal^  the  process  will  be 
reaction- rate  controlled  rather 
than  diffusion  controlled  when 
the  heating  rates^  and  conse¬ 
quently  the  surface  temperatures^  are  sufficiently  low.  If  mechanical  forces 
such  as  shear  stresses^  differential  pressures^  or  thermal  stresses  are  a  fac¬ 
tor  in  char  removal^  these  effects  must  be  superposed  on  the  effects  of 
oxidation. 

In  the  present  analysis  it  has  been  assumed  that  pyrolysis  occurs  at  the 
desired  back  surface  temperature  and  therefore  no  additional  insulation  or 
cooling  is  required.  Since  this  condition  is  not  usually  satisfied  with  actual 
materials^  the  analysis  can  be  used  only  to  calculate  the  amount  of  material 
that  will  be  pyrolyzed.  The  effect  of  the  heat  conducted  from  the  pyrolysis 
zone  into  the  uncharred  material  on  the  rate  of  pyrolysis  can  be  accounted  for 
by  adjusting  the  heat  of  pyrolysis  to  include  the  heat  required  to  raise  the 
material  from  its  initial  temperature  to  the  pyrolysis  temperature.  The  actual 
temperature  response  of  the  uncharred  material  and  of  the  structure  must  be 
determined  on  the  basis  of  a  separate  analysis. 

In  deriving  the  basic  differential  equation^  no  assumptions  have  been  made 
regarding  the  functional  dependence  of  the  enthalpy  or  of  the  radiant  heating; 
therefore^  if  these  parameters  are  expressed  as  functions  of  equation  (19) 

is  still  a  valid  relation^  and  solutions  can  be  obtained  by  nimmerical  methods. 


Figure  9.~  Heat  input  required  for  char  weight 
to  reach  l/e  of  extremum  char  weight. 
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In  particular^  if  the  enthalpy  and  radiant  heating  are  slowly  varying  func¬ 
tions  of  a  solution  can  be  obtained  by  successive  applications  of  equa¬ 

tions  (20a)  and  (22);  that  is^  the  final  value  of  ¥c  in  each  calculation 
becomes  the  initial  value  in  the  succeeding  calculation.  The  appropriate  cur¬ 
rent  values  of  all  the  environmental  parameters  must  be  used  in  each  step  of 
the  calculation. 


CONCLUDING  REMARKS 


A  diffusion- controlled  oxidation  mechanism  for  the  removal  of  char  at  the 
surface  of  a  char-forming  material  can  be  incorporated  in  an  approximate  anal¬ 
ysis,  and  the  performance  of  such  materials  in  various  thermal  environments 
can  be  evaluated.  For  a  given  material  and  environment,  there  exists  an 
extremum  char  weight.  The  performance  of  a  material  at  any  time  is  a  function 
of  the  ratio  of  the  actual  char  weight  to  the  extremim  char  weight. 

When  the  actual  char  weight  approaches  the  extremum  char  weight,  the  rate 
of  pyrolysis  approaches  a  limiting  value  which  is  independent  of  the  thermal 
conductivity  and  surface  emissivity  of  the  char  and  of  the  heat  of  pyrolysis. 
The  extremim  char  weight  increases  linearly  with  the  enthalpy  level;  therefore, 
the  actual  char  weight  approaches  the  extremum  value  much  more  quickly  in  tests 
conducted  at  low-enthalpy  levels  than  in  a  high-enthalpy  reentry  environment. 
This  factor  must  be  considered  in  evaluating  the  significance  of  results 
obtained  in  facilities  which  provide  a  low- enthalpy  test  stream. 

The  effects  of  radiant  heating  on  the  performance  of  charring  ablators 
appear  in  the  governing  equations  only  in  the  extremum  char  thickness  rela¬ 
tion.  The  analysis  indicates  that  with  constant  total  heat  input  the  effec¬ 
tiveness  of  such  materials  increases  with  the  substitution  of  a  small  component 
of  radiant  heating  for  convective  heating.  This  increase  in  effectiveness  is 
much  more  pronounced  at  low  enthalpies  than  at  high  enthalpies. 


Langley  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Langley  Station,  Hampton,  Va.,  December  5^  I963. 
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APPENDIX 


ANALYSIS  OF  OXIDATION 


The  oxidation  of  carbon  surfaces  in  hot  airstreams  has  been  extensively 
studied  analytically.  (See  refs.  5  to  7.)  It  is  found  that  at  low  surface 
temperatures  the  rate  of  oxidation  of  carbon  is  controlled  by  the  reaction 
rates.  At  higher  surface  temperatures,  however,  the  process  is  diffusion  con¬ 
trolled;  that  is,  the  rate  of  oxidation  is  limited  by  the  quantity  of  oxygen 
diffusing  to  the  surface  through  the  boundary  layer.  For  porous  carbon  sur¬ 
faces,  the  temperature  at  which  the  process  becomes  diffusion  controlled  is 
lower  than  it  is  for  nonporous  surfaces. 

The  char  layer  that  is  formed  when  charring  ablators  are  heated  tends  to 
be  porous,  as  evidenced  by  the  fact  that  the  gaseous  products  of  pyrolysis 
transpire  throiigh  it.  In  addition,  injection  of  these  pyrolysis  products  into 
the  boundary  layer  reduces  the  quantity  of  oxygen  diffusing  to  the  surface .  It 
is  anticipated,  therefore,  that  the  oxidation  of  the  char  surface  will  be  dif¬ 
fusion  controlled,  even  at  relatively  low  temperatures., 

The  rate  at  which  oxygen  diffuses  to  the  surface  is  determined  from  the 
boundary-layer  equation  for  conservation  of  oxygen.  This  equation  is  (ref.  8) 


p'u 


8x  8y  ay\Nsc 


The  energy  equation  can  be  expressed  in  the  following  form: 


(Al) 


1  \8u^ 
Npr/^ 


(A2) 


Equations  (Al)  and  (A2)  are  based  on  a  nonreacting  mixture  of  two  gases  having 
the  same  heat  capacity.  If  =  1  or  if  5u^/Sy  is  small  compared  with 

Sn/Sy^  which  is  valid  at  stagnation  regions^  equation  (A2)  may  be  expressed 
approximately  as  follows: 


P 


'uM+  p-v  M 
8x  oy 


5  /  H 

8y\Npj,  8y/ 


(A5) 


The  maximum  rate  of  diffusion  of  oxygen  to  the  surface  is  obtained  when 
the  concentration  of  oxygen  at  the  wall  vanishes.  The  boundary  conditions  on 
equations  (Al)  and  (A3)  are 


C  =  0 
H  = 


(A4) 


(y  =  0) 
(y  =  0) 
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c  c, 


H  H, 


(y  =  “) 
(y  =  «>) 


The  oxygen  concentration  and  the  enthalpy  can  be  replaced  by  the  dimen¬ 
sionless  variables: 

C  =  —  (A( 


ii  _  H  -  Hy 


(A6b) 


In  terms  of  these  new  variables,  equations  (Al),  (A3),  (a4),  and  (A5) 


p'u|^  +  P'v^  =  ^ 

dx  dy  dy\Nsc  dy 


(A7a) 


I  dH  ,  I  dH  _  d  /_jx  dH 
P  ^  dx  P  dy  dyUpr  dy 


C  =  H  =  0 


C  =  H  1 


(y  =  0) 
(y  =  “) 


(A7b) 


(A7c) 

(A7d) 


Because  of  the  formal  similarity  of  eq.uations  (A7a)  and  (A7b)  and  their 
boundary  conditions  (A7c)  and  (A7d)^  if  the  solution  of  equation  (A7t>)  is 


H  =  cp(x,y,Npp) 

the  solution  of  equation  (A7a)  is 

C  =  Cp(x,y,Nsc) 


(A8a) 


(A8b) 


where  cp  is  the  same  function  in  both  cases.  It  is  shown  in  reference  10 
that  cp(x^y^Npj.)  can  be  expressed  as 


cp(x,y,ITpr)  =  Np^Vi(x,y) 


(a8c) 


Therefore^ 


^(x^y.Ngc)  =  Ng- Vi(x,y) 


(A8d) 
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